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SUMMARY 


An analysis is presented to compute the film thickness, 
pressure and load relations between a rigid ball and rigid 
groove in normal approach when lubricated by a fluid with an 
exponential pressure-viscosity relationship. The geometry of 
the ball-groove system is reduced to the equivalent system of 
a paraboloid approaching a flat plate. 

Exact and approximate solutions are presented for the load 
and pressure relations. There is found to be a limiting load 
for a given geometry and lubricant regardless of the rate of 
approach. 


INTRODUCTION 


Comparatively little work has been done on squeeze films 
between curved surfaces. In reference 1, the squeeze film 
between a ball and spherical seat was investigated for the 
case of a constant viscosity lubricant. In reference 2 the 
elastohydrodynamic lubrication of shperes in normal approach 
was investigated using an exponential pressure-viscosity 
relationship. In reference 2, an analytical solution for the case 
of rigid surfaces and a numerical solution for deformable surfaces 
up to moderate pressures and for large film thicknesses are 
presented. For the range of pressures and film thicknesses 
considered in reference 2 it was found that the rigid surface 
solution yields results close to the elastic solution. 
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The work reported herein extends the analysis of reference 2 
to the case of a sphere approaching a cylindrical groove. The 
direction of approach is normal to the two surfaces. 


ANALYS IS 


The system under consideration together with the coordinate 
system used is shown in figure 1. In accordance with standard 
notation, the radius of curvature is considered positive when 
the center of curvature lies within the body. 

If the film thickness is much smaller than the radii of 
curvatures, then the customary parabolic approximation may be 
made and the system is approximately equivalent to a paraboloid 
approaching a flat plate as shown in figure 2. The equivalent 
radii of curvature in the x-z and y-z planes are given by 

R R„ 

R - GX 

x - R + Rp 
Gx 


R R 


R = 
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R + R 


2iL 


Gy 


where the signs of R Qx and R Gy are taken as defined previously. 
Using the parabolic representation, the film thickness at x y 
is given by 

2 2 


h 




= h Q + Ax 2 + By 2 


( 1 ) 
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The Reynolds equation for the lubricant film is given in 


reference 3 ° 



'Hie pressure-viscosity relationship will be considered 
to follow the exponential law: 


p 

A " A e 


(3) 


This relationship has been shown in reference 4 to be invalid 
at high shear when the pressure exceeds 3*8 10 N/m ( 55 000 psi). 

However, in reference 2 it is found that elastic considerations 
become important as the central pressure increases beyond a 
value equivalent to oC p Q - 5 

— 8 p 

For a synthetic paraffinic lubricant (c* = 1,3 10 m /N) 

8 2 

this corresponds to a central pressure of 3*86 10 N/m (56 000 psi). 

Thus, the breakdown of the exponential pressure-viscosity 
relationship and the rigid surface approximation occur at about 
the same pressure. 

Substituting the pressure-viscosity relationship (3) in 
the Reynold's equation ( 2 ) and considering the lubricant to 
be isothermal, the following partial differential equation 
is obtained 


a_ ( k 3 ip. 

axle " 13 <>*/ d y 


(4) 
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The particular solution to this equation is 


P 




(A + B)h 2 


letting 


Vo h « 

(A + B) 


= -k‘ 


this equation becomes 


p = 




( 5 ) 


(6) 


In accordance with the arguments presented in reference 5 
and reference 6 that the pressure distribution is only weakly 
dependent upon the shape of the boundaries, equation (6) will 
be considered as the complete solution. 

The normal load is obtained by integrating over the area 
and since the pressure drops off rapidly with distance from 
the origin, the upper limits on x and y will be taken as infinity. 
Therefore the load W is given byi 



( 7 ) 


Integrating (7), the expression for W is found to be 
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A simpler approximation for the pressure and load may be 
obtained by writing the pressure as a series i 



For small values of kA, higher order tens can be neglected 
and therefore i 


p = 



( 9 ) 


Integrating as before i 



0 O 


dy 


= _JL 
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( 10 ) 


The solutions may be written in dimensionless form as 
follows i 
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( 11 ) 
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and for the approximate expression i 
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A comparison of the exact and approximate solutions for the 
load parameter versus central film thickness is shown in figure 3* 
The limiting case for the exact expression as h Q -> k isi 


Woe 

2 rr h o j Ry R y 


In 4 = 1.386 


This is the same value obtained by Christenson (ref. 2) for 
the case of a sphere and plate. 

As in the case of the sphere on a flat plate this implies 
a limiting load for a given fluid and geometry which is given 
by the expression 


w = 2. 772 7 t/ R x R y K , 

OL 

The pressure distribution may also be written in dimensionless 
form by use of the variables: 


H = h/h Q 
H = h A 0 
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then 


K = kA 0 

P S «tp 
X = x//z R^h" 

y - y//2s/“ 

H = 1 + X 2 + Y 2 

and for the exact solution 

P = - In (l - K 2 /!! 2 ^ 
and farther approximate solution 

P = K 2 /H 2 

Typical pressure distribution along the major axis for the 
exact and approximate solutions for K = 0.7 are shown in 
figure 4. 

CONCLUDING REMARKS 

The pressure distribution follows the same trend as for rigid 
spheres in normal approach and differs considerably from the 
Hertzian pressure distribution. As in the case of spheres, 
there is a limiting load which is reached when the central 
pressure approaches infinity. In reality, of course, significant 
elastic deformation will occur before this load is attained and 
the pressure distribution will be altered. 
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APPENDIX - SYMBOLS 


A 

B 

H 

h 

h 

K 

k 

P 


P 

Po 

R 

R 


R 


Gx 

G y 


w 

x 

Y 


x,y,z, 

o' 


l/2R x , m _1 (in) 
l/2Ry, m”* (in) 

dimensionless film thickness * hAg 
film thickness, m (in) 
central film thickness, m (in) 

k Ap. . 

J - 3yu 6 ho^/ (ft + B), m (in) 

dimensionless pressure 
2 

pressure, N/m (psi) 

p 

central pressure, N/m (psi) 
radius of ball, m (in) 

radius of curvature of groove in x-z plane, m (in) 

radius of curvature of groove in y-z plane, m (in) 

radius of curvature of paraboloid in x-z plane, ra 

radius of curvature of paraboloid in y-z plane , m 

normal load, N (lb) 

dimensionless coordinate = x/Jz hgR x 

dimensionless coordinate = y//2h 0 Ry 

Cartesian coordinate, m (in) 

pressure-viscosity exponent m /N (psi -i ) 

2 2 

absolute viscosity, N-sec/m (lb-sec/in ) 

2 2 

ambient viscosity, N-sec/m (lb-sec/in ) 


(in) 

(in) 
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Fig. 3 COMPARISON OF APPROXIMATE AND EXACT SOLUTIONS 
FOR LOAD VARIABLE 
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Fig. 4 EXACT AND APPROXIMATE PRESSURE DISTRIBUTIONS 
FOR k = 0.7 
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